Dissociation of quarkonium in an anisotropic hot QCD medium 
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Abstract 

We have investigated the properties of quarkonium states in an anisotropic hot 
QCD medium by correcting the fuh Cornell potential, not the Coulomb term alone as 
usually done in the literature, with a dielectric function from the hard-loop resummed 
gluon propagator. We have found that in-medium modification in anisotropic medium 
causes less screening than in isotropic medium. In the short distance limit, potential 
does not show any medium dependence whereas in the long-distance limit, it reduces 
to the Coulomb potential with a dynamically screened color charge. In addition, 
anisotropy in momentum space introduces a characteristic angular dependence in the 
potential and as a result, quarkonium states in anisotropic medium are more tightly 
bound than in isotropic medium. In particular, quark pairs aligned in the direction 
of anisotropy are more bound than perpendicular to the direction of anisotropy. 
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1 Introduction 

Ultrarelativistic heavy-ion experiments have shown very rich physics which cannot be inter- 
preted by mere extrapolation from elementary nucleon-nucleon collisions to nucleus-nucleus 
collisions. This is evidenced from the suppression of high transverse momentum region of 
hadron spectra up to a factor of 5 relative to nucleon-nucleon collisions, which is an indi- 
cation for strong absorption of high-energy partons traversing the medium [1]. Also the 
inclusive production of charm quark bound states are suppressed by a factor of 3-5 at both 
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Super Proton Synchrotron (SPS) and Relativistic Heavy Ion Collider (RHIC) experiments 
which hints their dissolution in the medium [2, 3]. These accumulated informations indicate 
that a new form of matter has been produced in the heavy-ion experiments called as quark- 
gluon plasma (QGP). Among different experimental observations which may be served as 
the signals for the QGP formation, heavy quarkonium dissociation has been proposed long 
time ago as a clear probe of the QGP formation in collider experiments [4, 5]. First, in a 
pioneering work by Matsui and Satz [6] and later in a follow up quantitative calculation [7], 
it was shown that the suppression of J/ip^s could be explained by the (color) screening of 
the potential between a heavy quark and anti-quark by the deconfined light quarks and 
gluons. Theoretically, one can study the quarkonium states by using effective field theories 
(EFT). Since the mass of heavy quark, rriQ is much larger than the intrinsic scale in QCD, 
Aqcd, the heavy quark and anti-quark are expected to move slowly with a relative velocity 
V <^ 1 and results non-relativist ic version of QCD (NRQCD) [8, 9]. However, NRQCD 
does not fully exploit the smallness of v which further gives rise potential non-relativistic 
QCD (pNRQCD) [10, 11] by integrating out modes above both rriQ and m^v. 

The heavy quark pairs formed in relativistic nuclear collisions develop into the physical 
resonances and traverses the plasma and then the hadronic matter before decaying into a 
dilepton pair. Even before the resonance is formed it may be absorbed by the nucleons 
streaming past it [12]. By the time the resonance is formed, the screening of the color force 
in the plasma may inhibit the formation of bound states. The resonance(s) could also be 
dissociated either by hard gluons [13, 14, 15, 16, 17] or by comoving hot hadrons [18]. In 
order to disentangle the above mentioned effects [19], we must know how the properties of 
quarkonium states change in medium. The basic tools of phenomenological approach to 
study the properties of quarkonium states are potential models where relativistic effects can 
also be incorporated for excited states of charmonium. At zero temperature, the potential 
model has made great success. At finite temperature, the essence of the potential model 
in the context of deconfinement is to use a finite temperature extension of the potential at 
zero temperature. Quantitative understanding of the bound state properties needs the exact 
potential at finite temperature which, in principle, should be derived directly from QCD, 
just as the Cornell potential at zero temperature. Such derivation at finite temperature, 
however, is plagued by the existence of temperature-driven scales, T, gT, g'^T [20, 21]. 
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Due to these difficulties, the lattice-based potentials become popular. However, neither the 
free energy nor the internal energy can be directly used as the potential. In fact, what 
kind of screened potential should be used in the Schrodinger equation which describes 
well the bound states at finite temperature is still an open question. However, recently 
more involved calculations of quarkonium spectral functions and meson current correlators 
obtained from potential models have been performed which have been compared to first- 
principle QCD calculations performed numerically on lattices [22, 23, 24, 25, 26, 27, 28]. 
Regarding the criteria of dissociation, besides the binding energy for a particular potential, 
the decay width is another important quantity to determine the dissociation of the bound 
states. The calculation based on a real-valued potential model does not include the true 
width of a state. By performing an analytical continuation of the Euclidean Wilson loop 
to Minkowski space, the potential has an imaginary part due to Landau damping which 
clearly broadens the peak [21, 29]. 

In the RHIC or LHC era (small /is), it has been confirmed from the recent studies that 
the transition from nuclear matter to QGP is not a phase transition, rather a crossover [30]. 
The large distance property of the heavy quark interaction is important for our understand- 
ing of the bulk properties of the QCD plasma phase, e.g. the screening property of the 
quark gluon plasma [31], the equation of state [32, 33, 34] etc. In these studies, deviations 
from perturbative calculations and the ideal gas behavior are expected and were indeed 
found at temperatures larger than the deconfinement temperature. It is then reasonable 
to assume that the string-tension does not vanish abruptly at the deconfinement point, so 
one should study its effects on heavy quark potential even above Tc. This issue, usually 
overlooked in the literature where only a screened Coulomb potential was assumed above 
Tc and the linear/string term was assumed zero, was certainly worth of investigation. Re- 
cently [35] a heavy quark potential at finite temperature was derived by correcting the 
full Cornell potential (linear plus Coulomb), not its Coulomb part alone, with a dielectric 
function encoding the effects of the deconfined medium. This was found to have an addi- 
tional long range Coulomb term, in addition to the conventional Yukawa term. In the short 
distance limit, the potential is reduced to vacuum potential, i.e., does not see the medium. 
On the other hand, in the large distance limit (where the screening occurs), potential is 
reduced to a long-range Coulomb potential, apart from a constant term. Thereafter the 
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binding energies and dissociation temperatures of the ground and the first excited states 
of charmonium and bottomonium spectra [35, 36] were determined. The estimates with 
the Debye mass extracted from the lattice free energy are consistent with the finding of 
recent works based on potential models [37, 38, 39], whereas the Debye mass in the leading- 
order gives the identical results based on the lattice correlator studies and also with the 
potential model studies with stronger binding potential, i.e. internal energy [40, 41]. This 
finding ensues a basic question about the nature of dissociation of quarkonium in a hot 
QCD medium. 

However, this study was limited to an isotropic quark-gluon plasma. Actually, the par- 
ton system generated in the ultra-relativistic heavy-ion collisions cannot be homogeneous 
and isotropic because the longitudinal expansion at early times is much faster than the ra- 
dial expansion. At weak coupling, the longitudinal expansion causes the system to quickly 
become much colder in the longitudinal direction than in the transverse direction. As a 
result, an anisotropy in momentum space sets in and causes the parton system produced at 
the early stage of collision unstable with respect to the chromomagnetic plasma modes [42]. 
These instabilities facilitate to isotropize the system and hence quicken the process of its 
equilibration [43, 44, 45]. Motivated with this preamble on momentum anisotropy, we have 
made an attempt to investigate its effect on the heavy-quark potential and subsequently 
on the dissociation of quarkonia states. 

Our work is organized as follows. In section 2, we will discuss the medium modifications 
to a heavy quark potential both in isotropic and anisotropic medium. Recently one of us 
studied medium modification in isotropic medium [35], so first of all, we will start with 
their work in Sec. 2.1 and then extend the discussion to a medium which exhibits a local 
anisotropy in momentum space in Sec. 2. 2. To do that, we first obtain the self-energy ten- 
sor for an anisotropic medium, unlike isotropic medium, to obtain the hard-loop resummed 
gluon propagator. Thereafter the dielectric permittivity tensor be defined in terms of prop- 
agator and its Fourier transform at vanishing frequency gives the desired non-relativistic 
potential at finite temperature which is found to have a dependence not only on the relative 
separation of QQ pair but also on their relative orientation with respect to the direction of 
anisotropy and always deeper than in the isotropic medium. With this potential at finite 
temperature, we solve the Schrodinger equation to explore the properties of quarkonium 
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states in medium and determine their dissociation temperatures, respectively, in Section 
(s) 3.1 and 3.2, respectively. Finally, we conclude in section 4. 



2 Heavy-quark effective potential 
2.1 For isotropic medium = 0) 

Potential models are based on the assumption that the interaction between a heavy quark 
and its antiquark can be described by a potential. At T = 0, the hierarchy of well separated 
energy-scales: mq ^ rnqv ^ rnQv"^ , allows one to systematically integrate out the different 
scales and obtain the non-relativistic potential QCD (pNRQCD) where the Cornell poten- 
tial indeed shows up as the zeroth-order matching coefficient [10, 11]. Inspired by its success 
at zero temperature, the potential model has been applied at finite temperature, with the 
main assumption that medium effects can be accounted by a temperature-dependent po- 
tential. 

Recent lattice results [30] indicate the phase transition in full QCD appears as an 
crossover rather than a 'true' phase transition with the related singularities in thermody- 
namic observables. In light of the above findings, one cannot simply ignore the effects of 
string tension between the quark-antiquark pairs beyond T^. This is indeed a very impor- 
tant effect which needs to be incorporated while setting up the criterion for the dissociation. 
Recently [35, 36], this issue has successfully been addressed for the dissociation of quarko- 
nium in QGP by one of us and we closely follow their work in brief. Let us now start with 
a heavy quark potential (Cornell potential) at T=0 



where a and a are the phenomenological parameters. The former accounts for the effective 
coupling between a heavy quark and its antiquark and the latter gives the string coupling. 
The medium modification enters through the Fourier transform of heavy quark potential 
as 
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where the Fourier transform (FT) of the Cornell potential, V{k) is given by [35] 

■ a 4a 



Vit)^-^)-^ (3) 



and the dielectric permittivity, e{k) is given in terms of the static limit of the longitudinal 
part of the gluon self-energy [46]. In isotropic case, it can be decomposed into two compo- 
nents longitudinal (Hl) and transverse (Ht) where, in the static limit, are associated with 
electric and magnetic screening but transverse part IlT{0,k — 0, T) vanishes i.e., static 
magnetic fields are not screened. In perturbation theory, the quantity that enters in the 
Fourier transform of the potential at finite temperature is the static limit of the longitu- 
dinal gauge boson self-energy. The perturbative one- loop thermal contribution to 11^ was 
calculated long ago as [47] 



U,iO,k,T) = g^T^^^ + ^j^ml , (4) 

which defines the screening mass m^. Since the static limit of the self-energy is momen- 
tum independent, the pole of the inverse of dielectric permittivity is simply the gauge 
invariant Debye mass m^, so it leads to an exponential damping of the potential V{r) ~ 
exp(— m^r)/r. In particular, this form of II^, has the consequence that gluons screen the 
strong interaction, in contrast to the zero temperature case, over long distances. If one 
assumes non-perturbative effects such as the string tension which survive even above the 
deconfinement point [48] then the dependence of the dielectric function on the Debye mass 
may get modified. However, we incorporate the nonperturbative effects through the Fourier 
transform of linear string term. Finally, one can define a dielectric permittivity in one-loop 

After substituting it in the Fourier transform of Cornell potential (2) and then evaluating 
its inverse FT, one obtains the r-dependence of the medium modified potential [35, 36] as 

VirT) = [J^e^-rm 
^ ' ' J {2Tifl^ e{k) 

f 2a \ e"^ 2a 2a 
= a — H , (6) 

J r m^r 

with f = m^r. The constant terms are introduced to yield the correct limit of V{r,T) 
as T — 7- 0. Such terms could arise naturally from the basic computations of real time 
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static potential in hot QCD [29] and also from the real and imaginary time correlators 
in a thermal QCD medium [49]. The medium modified potential thus obtained has an 
additional long range Coulomb term with an (reduced) effective charge, in addition to the 
conventional Yukawa term. 

In the small distance limit r -C l/m^, the above potential reduces to the Cornell 
potential, i.e. no medium effects. On the other hand, in the screening region r ^ 
the potential (6) reduces to 

V{r) ~ — ^ - am^ (7) 

D 

The above form (apart from a constant term) is a Coulombic type as encountered in hydro- 
gen atom problem with identifying the fine structure constant with the effective charge 
2cr/m^. The binding energies and dissociation temperatures for quarkonium states can 
thus be determined by solving the Schrodinger equation numerically with the full poten- 
tial (6) or analytically with the approximated form (7). In addition, one can also exploit 
the advantage to demonstrate the fiavor dependence of the dissociation process where the 
dissociation temperatures for 2-fiavor are found to be higher than the 3-fiavor case [35, 36]. 

2.2 For anisotropic medium ^ 0) 
2.2.1 Dielectric permittivity tensor 

To study the perturbative potential with an anisotropic parton distribution, consider a hot 
QCD plasma which, due to expansion and finite (momentum) relaxation time, manifests a 
local anisotropy in momentum space through the distribution function 



i.e., f(k) is obtained from an isotropic distribution /jso(|k|) by removing particles with a 
large momentum component along the direction of anisotropy [42], n. We shall restrict 
ourselves to a plasma close to equilibrium and so that /iso(k) is either a Bose- Einstein 
distribution nsik) or a Fermi-Dirac distribution npik) function. This may be true because 
by the time quarkonia have been formed from the initially produced QQ pairs, system may 
not be highly anisotropic rather closer to isotropic distribution. So, in the limit of small 
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anisotropy, anisotropy parameter ^ is related to the shear viscosity by 

(9) 

ITS 

where 1/r denotes the expansion rate of the fluid element and 77/5 is the ratio of shear 
viscosity to entropy density. However, the degree of momentum-space anisotropy is also 
determined by the parameter, 

j]4)__ 

where ki = k.n and kj^ = k — n(k.n) are the components of momentum parallel and 
perpendicular to the direction of anisotropy, n. The positive value of anisotropic param- 
eter [C, > 0) corresponds to the squeezing of the distribution function in the direction of 
anisotropy and — 1 < < corresponds to the stretching in the direction of anisotropy. 
The calculation of the real part of the potential at finite anisotropy was first obtained in 
Ref. [29, 42, 50] and was later extended to include the imaginary part [21, 49, 51, 52, 53]. 
To study the effect of anisotropy on the medium modified potential we have to find the 
self energy in anisotropic medium. With the specified anisotropic distribution function, we 
can compute the gluon self-energy analytically [54]. We will restrict our consideration to 
the spatial part of U'^'^ for simplicity and the time-like components can be easily obtained 
by using the symmetry and transversality of the gluon self-energy tensor. The spatial 
components of the self-energy tensor reads [50] 

where = {po,p) is the four momentum of external gluon, p = |p| = amplitude of spatial 
momentum. (I7 v = k/|k|) is a light-like four vector with v = |v| and df /dk\ in terms 
of new variable /c, is given by 

dm_ ^ v' + i{^.ny dfCk') .^2) 
dk^ v/l + e(v.n)2 d~k ' 

where k'^ = ^^(1 + ^(v.n)^). After integrating out k, self energy, Uij simplifies into [50] 

^^^-"^-J in"" {l+a^r.nyyy ^ P■V + ^e) ' ^^^^ 
where the square of Debye mass, m? is defined by 



00 



= / dkk'^^^^. (14) 



Unlike isotropic medium, in anisotropic medium n'j shows an extra dependence on the 
preferred anisotropic direction (n), therefore it can no longer be decomposed into transverse 
and longitudinal parts rather it becomes a tensor [42, 55, 56]. So the self-energy tensor be 
decomposed into four structure functions as [42] 

n*^' = aA'^ + (5B'^ + -iC'^ + 5D'^ , (15) 

where the coefficients a, /3, 7 and 5 can be determined for any value of anisotropy parameter 
^ [42]. Using this, one can find the gluon propagator for an anisotropic medium in temporal 
axial gauge as [42] 



(A;2 - w2 + + 7)(/3 - a;2) - ^2^252 

In order to see how the anisotropy affects the response to static electric field, we examine 
the propagator in the static limit (a; — ?■ 0). Defining the masses in the static limit [42] 

mi = lim a, mi = lim -13 

° uj->-o P a;2 

mi = lim 7, m? = lim Im S , (17) 

the gluon propagator becomes 

+ ml + mi) kikj 



Aij{uj, k) 



Cj2 [(/^2 _|_ ^2 _|_ jy^T^(J^2 _|_ j^2-j _ j^4j 

We can now factorize the denominator of the gluon propagator as 

+ ml + m^) {e + mj) - = ( + ml) (k^ + mi) (19) 

where 



2ml = ± y^M^ - 4(m2 (m2 + m2) - m^) , = ml + m^ + m"^^ (20) 

Thus the dielectric permittivity for anisotropic medium in temporal axial gauge can now 
be obtained from the definition [57] 

-i/M r 2hkj. , k\e + ml + m^) 

' (^) = -a" (P + o(P + m^_) ' ^21) 
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If the anisotropy is small, the pole masses (17) simplifies, in terms of anisotropy parameter 
into [42] 



ml = -|(1 + cos2/3„) m 
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(^l + ^(3cos2/3„-l)^ ml 



2 ^-20 2 

m^ = - sm (3n m^ 

m] = sin I3n cos I3n ml (22) 

where /3„ is the angle between k and n . So the explicit dependence of m± on anisotropy, 
in the small ^ limit is given by 

2 / 1 I ^ „„„ o/Q 1 N 1 ™2 



= (^1 + ^(3 cos2/3„-l)Jm 
= -| cos2/3„m2. (23) 

In the isotropic limit, all masses become zero except one (m^ = = = = 0, = 
ml), which is the only pole in the isotropic medium (5). 

2.2.2 Medium- modification to heavy quark potential 

Once we have obtained the dielectric permittivity in anisotropic medium (21), we substitute 
it in the Fourier transform (2) and then evaluate its inverse Fourier transform to obtain 
the medium modified potential [58] in an anisotropic medium: 



1 f ,3,. V{k) 



(27r)^/^ J e{k) 



a 



(P + m^ + m^) 



27r2 J (^k^ + ml){k^ + ml) 

^d^k „'f+';°;'"'' (24) 



(27r)2 J k\k^ + ml){k^ + mlj 
After substituting m+ and m_ from (23), the potential (24) becomes 

VUeT) =-— f e"- - 

271^ J k^ + ml{l + l{3cos2f3n-l)) 
4a r d^k e^"" 



(27r)2y P [P + m2 (l + |(3cos2/3„-l))] 
Vi(r,e,r) + r2(r,e,r) , (25) 
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where Vi{r, ^, T) and V2(r, ^, T) are the medium- modified potential corresponding to Coulom- 
bic ( — 1/r) and hnear/string (err) term, respectively. Both terms can be separated into 
isotropic and anisotropic contributions as 

V,{r,^,T) = vl'\r,^ = 0,T) + v}'\r,^,T) (26) 
V2{r,^,T) = vi'\r,^ = 0,T) + vi'\r,^,T) (27) 

where Vi^\r, ^ = 0, T) and Vi'^\r, ^, T) are the isotropic and anisotropic contributions, due 
to medium modification of the Coulomb term, respectively. Similarly, ^^^^(r, ^ = 0,T) and 
V2 (r, ^, T) are the corresponding terms due to to medium modification of the linear term, 
respectively. Let us now calculate them one-by-one. 

The isotropic part, Vi^\r,^ = 0,T) of the Coulomb term (along-with a constant term 
which was earlier introduced to yield the correct limit at T — )■ as discussed in Eq.(6)) is 
given by 

V}'\r,^ = 0,T) = -^e-'-am, , (28) 

and the anisotropic part, Vi'^\r,^,T) is given by 

yi'i^.i.T) . -i ^ ,nl j (? - cos^.„) . (20, 

One immediately observes that in anisotropic medium, in addition to inter-particle distance, 
an angular dependence appears in the potential. There are two cases which are especially 
interesting for the discussion. First, we consider the case where r is parallel to the direction 
of anisotropy, n. In this case Vi^\y,^,T) reduces to 

00 7r 

00 ° 

= e^/^^(r,e = 0,T)(2^-^-^-l) (30) 

where we have taken the anisotropy n along 2;-axis and the angle between 2;-axis and k is 9. 
Thus, the complete in-medium modification to the Coulomb term, for parallel alignment 
becomes 

Vi(r||n,e,r) = Vl^\r,i = Q,T) + V^^\v\\^,i,T) 

am„ - am„ - f — 1 2 f \ 

= r^e-^-am^-e— r^e-'- 2— 31 

r r \ r b J 
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Next we consider the other case when r is transverse to the direction of anisotropy then 
in-medium-modification to the Coulomb term, V-['^\t,^,T) becomes 

OO TV 2-77 



y/')(r±n,e,T)=-e 



a 



27r 



^r/')(r,e = o,r) 



-m 
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k'^dk sine dO 
(A;2 + m|)2 



Akr cos d 



COS 



^'n) sin^ e 



\-e^ 1 1 f 

— + T + 77 + 77 



(32) 



r 2 3^ 

where we take r along 2-axis and n lying in x|/-plane, so is the azimuthal angle and 0„ 
is the angle between n with x-axis. Thus, the complete in-medium modification to the 
Coulomb term, for transverse alignment becomes 



Vl(r±n,e,T) = \//')(r,e = 0,T) + y,^^^(r±n,e,T) 

'e^ - 1 1 



— - — e — am^ + 4 — - — e 



f 1 
^ ~ 3 ~ 2 



(33) 



Next we calculate the medium- modification to the linear term, V2(r,,^,T). The isotropic 
part as already calculated in Eq.(6) is given by 



^('Hr,e = 0,T) 



2a 



2a 2a 



(34) 



whereas the anisotropic contribution for r parallel to the direction of anisotropy is given by 
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Thus, the complete in-medium modification to the linear term, for r ||n becomes 



(35) 



y2(r||n,e,r) = Vi'\r,^ = 0,T) + V,''>{r\\n,^,T) 

2a 



2a 2a 
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However, when r is transverse to the direction of anisotropy, n, the anisotropic contribution 
becomes 



^(')(r±n,e,r) = -e 
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Thus, the complete in-medium modification to the hnear term for r _L n becomes 
V2{r±n,^,T) = vi'\r,^ = 0,T) + vl'\r±n,^,T) 



- e 
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(38) 



Therefore, the full medium-modified potential consisting of linear and Coulombic term 
Eqs. (31 and 36), for the quark pairs aligned to the direction of anisotropy, yields as 



^(r||n,^,T) = { — -am^)^ 



2a 



Aa 



2a 

m 



am. 



D 

2 2 



r 



r 

12 



am 



2 

f 



r 
6 



1 



The potential thus obtained has a long range Coulombic tail in addition to the standard 
Yukawa term. In the short distance limit (r ^ V'^d '^^ rm^^ <^ 1), it reduces to the 
vacuum potential (Cornell), i.e. QQ pairs are not affected by the medium. On the other 
hand, in the long-distance limit (r ^ 1/''^d rm^^ ^ 1) (where the screening occurs), we 
can neglect the Yukawa term and for large values of temperatures, the product am^ will 
be much greater than 2a /m^. Thus the potential (39) finally simplifies to the following 
form 



^(r II n,e,T) 



2a 
m^ r 



^ I — am^ 



(40) 



which shows that potential for QQ pairs aligned in the direction of anisotropy is stronger 
than in isotropic medium (7). 

On the other hand, when the quark pairs are aligned transverse to the direction of mo- 
mentum anisotropy (r ± n), the medium modification to the linear and Coulombic terms 
together (Eqs. (38) and (33), respectively) gives rise to the following form: 

e"^ 2a 2a 
— H am 

r m^r m^ 

4cT _^f{e^-l) (e^-7) 1 f 

f 6 



nr±n,e,T) = (— 
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In the short distance hmit, the potential is similarly reduced to the vacuum potential 
whereas in the long-distance limit (rm^ ^ 1), the form (41) is simplified to the form 

V{v ± n, T) = f 1 + i V «^z,, (42) 

which again shows that potential even for transverse alignment is still stronger for nonzero 
values of anisotropy than in isotropic medium. 

We have noticed overall that in the short distance limit, potential have not been af- 
fected by the anisotropic medium too. On the contrary, unlike the the potential in isotropic 
medium, anisotropy introduces an angular dependence in the potential and gives rise a char- 
acteristic dependence between the relative separation (r) and the direction of anisotropy 
(n). This is evidenced from the observation that potential in anisotropic medium is always 
stronger than in isotropic medium. In addition, quark pairs aligned along the direction 
of anisotropy feel more attraction than the transverse alignment because the inter-quark 
potential along the direction of anisotropy is less screened than transverse alignment. 

To see the effects of anisotropy on the heavy quark effective potentials (39) and (41), 
we have plotted Figures 1(a) and 2(a), for quark pairs (r) aligned along and perpendicular 
to the direction of anisotropy, n, respectively. First observation common to all figures is 
that interquark potential in anisotropic medium is always more attractive than the isotropic 
medium. This can be understood as follows: customarily, in the small anisotropic limit, the 
anisotropic distribution function may be obtained from an isotropic distribution /jso(|k|) 
by removing particles with a large momentum component along the direction of anisotropy 
(n) i.e. fiso{^J^^ + ^(k.n)2. This transpires in the reduction of number of partons in 
anisotropic medium than in isotropic medium i.e. naniso(0 = '^iso/ + i- Therefore, the 
(effective) Debye mass in anisotropic medium is always less than in isotropic medium and 
hence potential will be screened less i. e stronger than in isotropic medium. 

The second observation is that the quark pairs aligned along the direction of anisotropy 
are stronger than aligned perpendicular to the direction of anisotropy. This can also be 
understood in terms of explanation of first observation. For the parallel alignment, compo- 
nent of momentum to be removed is higher than transverse alignment so the distribution 
function for the former case in the higher momentum region will be reduced compared to 
the latter case. Therefore the Debye mass for parallel alignment will be slightly smaller 
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Figure 1 ; Medium- modified heavy-quark potential as a function of f (= rm^ ) for quark pairs parallel to 
the direction of anisotropy, n. Left panel (a) represents the full potential (39) divided by {g^CFm^) and 
right panel (b) represents the contribution of coulomb and string part in the full potential as a function of 
r. 
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Figure 2: The notations are same as in Figure 1(a) and 1(b) but for transverse alignment. 
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than transverse alignment (see the definition of Debye mass (14)) and hence potential will 
be less screened compared to transverse case. However the difference will not be much 
because the contributions in the Debye mass from partons having higher momenta are very 
small. 

To understand the effect of medium modification of the linear (string) term on medium 
modified heavy quark potential quantitatively, we have plotted separately the medium 
modification to the linear term, Coulomb term and the full Cornell potential in Figure 1(b) 
and 2(b), for parallel and transverse case, respectively. Medium modification to the Cornell 
potential contains two parts: one is due to medium modifications of linear term (err) and 
the another one is due to the medium modifications of attractive Coulomb term {-1/r). 
As usually done in the literature, medium modification to the linear term does not arise 
because string tension was assumed to be zero [26, 59, 60, 61, 62, 63] beyond deconfinement 
temperature [64]. Since string tension is found to be nonzero at Tc rather it approaches zero 
much beyond Tc and hence the medium modification to the linear term may be nonzero 
even at temperatures beyond Tc, although it is very small. In isotropic case, medium 
modification to the linear term remains positive up to 2-3 T^, making the potential less 
attractive compared to T = 0. On contrast, in anisotropic case medium modification 
to the linear term becomes negative and the overall full potential (including the medium 
modification of -1/r term) becomes more attractive compared to isotropic case. 

3 Properties of Quarkonium in medium 
3.1 Binding energy 

Apart from the constant term, the simplified form of the potential for heavy quark pairs 
aligned along and transverse in the direction of anisotropy (n) (Eqs. 40 and 42) resem- 
bles with the potential in the hydrogen atom problem, after identifying 2(7/m^ with the 
fine structure constant e^. Thus, the solutions of the Schrodinger equation for both the 
potentials are known to us and gives the energy eigenvalues for the ground states and the 
first excited states for charmonium {J/ip, ip' etc.) and bottomonium (T, T' etc.) spectra. 
Even if we take the full form of the potentials (39 and 41), Schrodinger equation can be 
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solved numerically. However, the results do not differ much with the simplified form of the 
potential 



Therefore, Schrodinger equation with potential (40) for QQ pair aligned along n gives 
the energy eigen values of bound states as 

4 = -|-«™„; £! = !^(n-l?)' . (43) 

where mq is the mass of the heavy quark. Apart from the constant term (— am^), — 
— -E'l/4, —Ej /9 etc. are the energies of QQ bound state, in the first, second, third etc. Bohr 
orbit, respectively. These energies are known as the ionization potentials or the binding 
energies for the nth QQ bound states which become temperature-dependent through the 
Debye mass and decrease with temperature. These binding energies are seen larger than 
in isotropic medium by an amount (1 + 4^/3)^. 

Similarly, for QQ pair aligned perpendicular to the direction of anisotropy, n, the solu- 
tion of the Schrodinger equation with the potential (42) is given by 

where, apart from —am^, —E^j-, —E^j-fA, —Ej-/9 etc. are the binding energies in first, 
second, third etc. Bohr state, respectively, for QQ pairs aligned transverse to n. These 
binding energies are also seen larger than in isotropic medium by an amount (1 + ^/6)^. 

Combining both results on the binding energies in anisotropic medium, quarkonium 
states are seen always bound more strongly than in isotropic medium. In particular, the 
binding is more for QQ pairs aligned along n than perpendicular to n. 

In our calculation, we consider the Debye mass obtained by fitting the (color-singlet) 
free energy in lattice QCD [38] (m'" ) as 



ml = lAmf, mf=9{T)T^^ + ^ (45) 

where m^^ is leading-order Debye mass where we use both one-loop and two-loop expression 
for QCD coupling constant, g(T) at finite temperature [62, 65, 66], to see the effects of 
running coupling on the dissociation process. 
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Figure 3: Variation of J/i(j binding energy (in GeV) with temperature for QQ pairs aligned along(left 
panel) and transverse (right panel) to the direction of anisotropy (n). 
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18 



Parallel alignment 



50 



4U I I I I I I I I I I I I I I I I I I I I I I I 
-\ 
- \ 

\ 

\ 



40 - ^. 



30 - 



20 - 
















... i= 


0.5 




^< 1= 


1.5 






0.9 



10 



o 



5 

-\ 

\ 

\ 

- \ 



6 ^ 



4 - 



2 - 



I'ransverse alignment 
I I I I I I I I I I I I I I I I I I I I I I I I 



— • 5= 





... 5= 


0.5 


5= 


1.5 




0.5 0.6 



0.7 0.8 0.9 
T/Tr 



Figure 5: Variation of T binding energy (in GeV) with temperature for QQ pairs aligned along(left panel) 
and transverse (right panel) to the direction of anisotropy (n). 
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Figures (3)-(4) show the variations of binding energies (in GeV) with temperature (in 
units of critical temperature) for the ground and first excited cc states. In each figure we 
have calculated binding energies for different values of anisotropy parameter ^, to see the 
effect of anisotropy on the binding energies, compared to the isotropic case. In the left (a) 
and right panel (b) of each figure, we have plotted binding energies for QQ pairs aligned 
along and perpendicular to the direction of anisotropy n, respectively. There are mainly 
three crucial observations common to all figures. First, as the anisotropy increases, the 
binding of QQ pairs get stronger with respect to their isotropic counterpart because the 
potential becomes deeper with the increase of anisotropy. This is due to the fact that the 
(effective) Debye mass in anisotropic medium is always smaller than in isotropic medium. 
As a result the screening of the Coulomb and string contributions is less accentuated and 
thus quarkonium states are bound more strongly than in isotropic medium. However, 
effects of anisotropy on the binding energies are not so pronounced on the excited states, 
compared to the ground states. 

Secondly, in all figures, there is a strong decreasing trend in binding energies with 
the temperature because screening becomes stronger with the increase of temperature, so 
the potential becomes weaker compared to T = and results in early dissolution of QQ 
pairs in the medium. Thirdly, the binding of QQ pairs which are aligned in the direction 
of anisotropy are seen much stronger than those aligned perpendicular to n because the 
interquark potential for the former case is more attractive than the later. Similar variations 
for bottomonia states (T, T') are shown in Figs.(5)-(6). Our results on the temperature 
dependence of the binding energies show a qualitative agreement with the similar variations 
shown in [67]. 

Thus a study of the temperature dependence of the binding energy is poised to provide 
a wealth of information about the dissociation pattern of quarkonium states in anisotropic 
thermal medium that can be used to determine the dissociation temperatures of different 
states in the next Section. 
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3.2 Dissociation temperatures for heavy quarkonia 

Dissociation of a two-body bound state in a thermal medium can be understood in the 
following way: when the binding energy of a resonance state drops below the mean thermal 
energy of a parton, the state becomes feebly bound. The thermal fluctuations then can 
destroy it by transferring energy and exciting the quark-antiquark pair into its continuum. 
The spectral function technique in potential models defines the dissociation temperature 
as the temperature above which the quarkonium spectral function shows no resonance- 
like structures but the widths shown in spectral functions from current potential model 
calculations are not physical. The broadening of states with the increase in temperature 
is not included in any of these models. In Ref.[38], the authors argued that one need not 
to reach binding energy (-Ebin) to be zero for the dissociation rather a weaker condition 
Eun < T causes a state weakly coupled. In fact, when E'bin — T, the resonance states 
have been broadened due to direct thermal activation, so the dissociation of the bound 
states may be expected to occur roughly around i^bin — T. Therefore the upper bound on 
the dissociation temperature (To) for the n-th bound state of QQ pairs aligned along the 
direction of anisotropy n is obtained from the relation 

\El\ = ^ + am^=Tn , (46) 

where string tension (cr) is taken as 0.184GeV^, and critical temperatures (Tc) is taken 192 
MeV for three-flavor QCD medium. It is found that the dissociation temperatures (T^'s) 
for quarkonia states are always minimum for isotropic case and increases appreciably with 
the increase of anisotropy > 0), viz. J/ip and ip' are dissociated at 1.61Tc and 1.34 Tc, 
respectively in isotropic medium = 0) while in anisotropic medium they are dissociated 
at much higher temperatures 2.35Tc and 1.83 T^, respectively at ^ = 1.5. Similarly the 
dissociation temperatures of bb states goes higher in anisotropic medium viz. at ^ = 1.5 the 
Td's of T and T' are 2.91 Tc and 2.23 Tc, respectively, corresponding to 1.92 Tc and 1.54 Tc, 
respectively in isotropic medium. This increase in dissociation temperatures is due to the 
fact that screening gets reduced in anisotropic medium because a fraction of distribution 
function with higher momentum does not contribute to the Debye mass. This, in turn leads 
to the enhancement of effective charge compared to the isotropic medium by an amount 
(1 + 4,^/3)^. Hence the potential becomes deeper and consequently binding becomes more 
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tight. The dissociation temperature for other values of anisotropy parameters are given in 
Table 1. 

Similarly the upper bound on the dissociation temperature (To) for the n-th bound state 
of QQ pairs which are aligned perpendicular to the direction of anisotropy, n is obtained 
from the relation 

\E^\ = ^ + am^=Tn , (47) 

where the dissociation temperatures for the lowest lying and first excited states of QQ 
pairs for transverse alignment are given Table 2. Like parallel alignment scenario, the 
dissociation temperatures increases slowly with the increase of anisotropy and are minimum 
for isotropic case. For example, J/ip's {ip"s) are now dissociated at higher temperature 
1.72Tc (1.41 Tc) at ^ = 1.5 compared to their isotropic result - 1.61Tc (1.34Tc) whereas 
T's (T"s) are dissociated at higher temperature 2.08Tc(1.65Tc) compared to the isotropic 
result. The excited states also survive at higher temperatures (Table 2) compared to 
isotropic results. These observations in transverse alignment can again be understood in 
terms of the enhancement of effective charge compared to the isotropic case by an amount 
(1 + ^/6). As a consequence, potential again becomes stronger than isotropic case and 
binding becomes more tight. However, in this case, increase is smaller than the earlier case 
(parallel alignment). 

So far we have studied the dissociation process with a temperature dependent QCD 
running coupling in one-loop. To see the effects of QCD running coupling, we have cal- 
culated again the dissociation temperatures with QCD coupling with two-loop which is 
smaller than one-loop value and the results are given in Table 3. The values of dissociation 
temperatures (Table 3) are seen slightly smaller compared to one- loop results (Table 1), 
at least, for lower values of anisotropy. This can be understood in the following way. The 
binding energy (e^/i/2) of the ground state in Coulomb like potentials (Bohr's theory) is 
proportional to the square of the fine structure (coupling) constant (e^) and the reduced 
mass (yu) (in this case it is mQ/2). The QCD coupling upto two- loop is slightly smaller 
than one-loop value so the binding energy in two-loop coupling is slightly smaller than 
one-loop interaction. So less thermal energy is needed to excite into the continuum and 
results in smaller values of dissociation temperatures. For example. To for J/il) is 1.61 Tc 
with the one-loop QCD coupling whereas with the two-loop coupling, becomes slightly 
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State 


e = o.o 


e = o.5 


^ = 1.0 


e=1.5 


J/^ 


1.61 


1.90 


2.14 


2.35 


ip' 


1.34 


1.52 


1.69 


1.83 


T 


1.92 


2.32 


2.64 


2.91 




1.54 


1.81 


2.03 


2.23 



Table 1: Dissociation temperatures (Td) for the quarkonium states in anisotropic medium with strong 
couphng upto one-loop. When quark pairs are ahgned along the direction of anisotropy(r||n) and Debye 
mass from lattice parameterized free energy {m^ = 1.4m^'-') is taken with Nf = 3. 



State 


e = o.o 


e = o.5 


^ = 1.0 


e=1.5 




1.61 


1.65 


1.69 


1.72 


^' 


1.34 


1.36 


1.39 


1.41 


T 


1.92 


1.98 


2.03 


2.08 




1.54 


1.58 


1.61 


1.65 



Table 2: Dissociation temperatures {Tjj) for the quarkonium states in anisotropic medium with strong 
coupling upto one-loop. When quark pairs are aligned transverse to the direction of anisotropy(r _L n) and 
Debye mass extracted from lattice free energy (m^ = 1.4to^'-') is taken with Nf ~ 3. 

smaller, i.e. 1.51 T^.. However, for T state at higher values of anisotropy = 1.5), the 
decrease in binding energy due to running coupling from one-loop to two-loop is not only 
compensated by the enhancement (1 + 4,^/3)^ of the effective charge and the reduced mass 
but may even surpass than the one- loop value. As a result, the dissociation temperature 
for T with the two-loop coupling (2.96 T^) is slightly higher than one-loop value (2.92 T^) 
which is opposite to smaller anisotropy result. However, for quark pairs aligned transverse 
to the direction of anisotropy, the dissociation temperatures (shown in Table 4) are always 
smaller while going from one-loop to two-loop because the compensation or enhancement 
factor (1 +^/6)^ is much smaller than earlier case. As a result, effect of QCD couphng 
from one-loop to two-loop is not much pronounced. 



State 


e = o.o 


e = 0.5 


^ = 1.0 


e=1.5 




1.50 


1.85 


2.12 


2.35 


^' 


1.16 


1.40 


1.60 


1.77 


T 


1.87 


2.31 


2.66 


2.96 




1.42 


1.74 


2.00 


2.22 



Table 3: Dissociation temperatures {To) for the quarkonium states in anisotropic medium with strong 
coupling upto two-loop. When quark pairs are aligned along the direction of anisotropy(r||n) and Debye 
mass from lattice parametrized free energy (m^ = 1.4771^*^) is taken with Nf = 3. 
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State 


e = o.o 


e = o.5 


^ = 1.0 


e=1.5 




1.50 


1.55 


1.60 


1.64 


ip' 


1.16 


1.19 


1.22 


1.25 


T 


1.87 


1.94 


2.00 


2.05 




1.42 


1.47 


1.51 


1.55 



Table 4: Dissociation temperatures (Tu) for the quarkonium states in anisotropic medium with strong 
couphng upto two-loop. When quark pairs are ahgned transverse to the direction of anisotropy(r _L n) and 
Debyc mass from lattice parametrized free energy (m^ ~ 1.47ti^'^) is taken with Nf = 3. 

In the above mentioned results shown in Table 1-4, we have used the lattice form of 
the Debye mass (45). However, if we use the leading-order Debye mass which is 1.4 times 
smaller than the lattice Debye mass, the screening of the potential will be much smaller. 
This increases the binding energies (1/m^) substantially and results in large increase of 
the dissociation temperatures. For example, for the parallel alignment, J/ip will now be 
dissociated at much higher temperature 2.85 Tc at ^ = 1.5 compared to the value (2.35 Tc) 
whereas T is also dissociated at much higher temperature 3.60 T^. However, for transverse 
alignment, the increase of dissociation temperature due to perturbative result of Debye 
mass is not so large because the screening factor (1 + ^/Q)"^ is smaller than the former. 



4 Conclusions and Outlook 



In conclusion, we have studied the dissociation phenomena of quarkonia by correcting 
the full Cornell potential with a dielectric function from the hard-loop resummed gluon 
propagator embodying the effects of anisotropic deconfined medium and not its Coulomb 
part alone as usually done in the literature. We have noticed that in-medium modification 
in anisotropic medium causes always less screening than in isotropic medium and as a result, 
potential gets stronger always. In the screening region, it reduces to the Coulomb potential 
with a dynamically screened color charge. Further, anisotropy introduces a characteristic 
angular dependence in the potential and as a consequence, potential for quark pairs aligned 
in the direction of anisotropy are more attractive than pairs perpendicular to the direction 
of anisotropy. 

We have then systematically studied the temperature dependence of binding energies 
for the ground and first excited states of charmonium and bottomonium spectra in a three 
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flavor QCD medium. We have found that quarkonia states in anisotropic medium are always 
more bound than in isotropic medium. As a consequence, quarkonia states survives much 
higher temperature, compared to isotropic medium. This can be understood quahtatively 
by the fact that as the anisotropy of the medium increases, the screening becomes weaker 
and hence binding gets stronger. In particular, QQ pair aligned along the direction of 
anisotropy is more bound than aligned transverse to the direction of anisotropy. 

Finally we have determined the dissociation temperatures which are are consistent with 
the finding of recent theoretical works based on heavy quark potential in an anisotropic 
plasma [50]. However, our results on dissociation temperatures are in the higher side due 
to medium modification of the linear term which was absent in their calculation [50]. In 
fact, one-dimensional Fourier transform of the Cornell potential in the medium yields the 
similar form used in the lattice QCD to study the quarkonium properties which assumes the 
one-dimensional color fiux tube structure [68]. However, at finite temperature that may not 
be the case since the fiux tube structure may expand in more dimensions [61]. Therefore, 
it would be better to consider the three-dimensional form of the medium modified Cornell 
potential which has been done exactly in the present work. 

In brief, J/ip is found to be dissociated at much higher temperature 2.35Tc at ^ = 1.5 
compared to the lower temperature 1.61Tc in isotropic medium [E, = 0). The corresponding 
dissociation temperatures for the ground state (T) of bb state are 2.91 Tc and 1.92Tc in 
anisotropic and isotropic medium, respectively. We have also studied the sensitivity of QCD 
running coupling on the dissociation process by going from one-loop to two-loop calculation. 
For smaller anisotropy, the dissociation temperatures are found to be smaller whereas for 
higher anisotropy, T survives higher temperature. However, when the perturbative result 
of Debye mass is used, the dissociation temperatures are found to be much higher because 
the screening of the potential is less due to smaller value of the Debye mass. This finding 
envisages a basic question about the nature of dissociation of quarkonium in an anisotropic 
hot QCD medium. 

Apart from the uncertainty of the correct form of the potential, there is an arbitrariness 
in the definition of dissociation temperature. So, in future, we wish to investigate the 
dissociation through the decay width of quarkonium bound states calculated from the 
imaginary part of the potential. 
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